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Abstract 

A generalized propagation matrix method is used to study how scattering off local Einstein 
phonons affects resonant electron transmission through quantum wells. In particular, the parity 
and the number of the phonon mediated satellite resonances are found to depend on the available 
scattering channels. For a large number of phonon channels, the formation of low-energy impurity 
bands is observed. Furthermore, an effective theory is developed which accurately describes the 
phonon generated sidebands for sufficiently small electron-phonon coupling. Finally, the current- 
voltage characteristics caused by phonon assisted transmission satellites are discussed for a specific 
double barrier geometry. 

PACS numbers: 73.22.-f, 73.22. Gk,61.46.+w,74.78.Na 
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I. INTRODUCTION 



Resonant tunneling through quantum wells has been extensively studied in semiconduc- 
tor heterostructures, such as GaAs/Al x Gax_ x As double barriers ]l[ 2], 3], ^,2]. More re- 
cently, analogous electron transmission processes have also been investigated in the context 
of molecular junctions , llOj and mesoscopic rings [ill . Il2j | . Resonant tunneling is 



a purely quantum effect whereby electrons pass through structures made of potential wells 
and barriers with unit or near-unit transmission probabilities if they enter the quantum well 
at the particular energies of the structure's bound states. Following the initial experimental 



observation of satellitepeaks of t 



ical work 
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rese transmission resonances [l| , a large volume of theoret 



25, |26J has focused on the effects of 



phonon scattering on the electronic tunneling. Early on, it was recognized that perturbative 
treatments tend to miss the essential feedback effects between elastic and inelastic channels 
which lead to these satellite features in the electronic transmission ljj. In particular, it was 
found that electron-phonon scattering processes can cause the formation of polaronic bound 
states, leading to phonon assisted resonant tunneling 0, H, Q, \v\ . More recent works 
have shown that within a tight-binding description these features are further enhanced 18) , 



and phonon bands can form 



22 



231 ] . Furthermore, theoretical models have been general- 



ized to include the effects of the three-dimensional environment 
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211 ]. non-equilibrium 



dynamics 20j, and finite temperatures 24]. 

In this paper, we examine the hierarchy of polaronic resonances in the electron trans- 
mission through quantum well structures. In particular, we focus on even-odd effects with 
respect to the number of available phonon channels and on the emergence of impurity bands 
as this number becomes large. We also apply an effective theory which reproduces the de- 
pendence of the resonance peaks on the electron-phonon coupling strength and the phonon 
energy in the limit of sufficiently small coupling. The method we are using is a generalization 
of the propagation matrix technique 23] which takes into account elastic electron scattering 
at potential steps as well as scattering off local Einstein phonons. This approach allows 
a numerically exact calculation of the electron transmission through quasi-one- dimensional 
heterostructures without any perturbative constraints, such as limitations to particular pa- 
rameter regimes, or restrictions to specific energy ranges, such as low-energy resonant states. 
In addition, our method accounts for the feedback between an adjustable number of phonons 
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and the elastic transmission channel, and is therefore suitable to accurately describe the in- 
terplay between non-perturbative resonances of the many-body system. 

Before proceeding to the discussion of resonant tunneling through specific semiconduc- 
tor double wells, let us briefly point out some similarities and differences of this system 
with electronic transport through molecular junctions. 0, Q, 0, 0, [h]] In the theory of both 
physical systems, many-body methods are combined with scattering theory to obtain the 
tunneling density of states and the resulting current-voltage characteristics for electronic 
transport through small objects with quantized energy levels. In both cases one observes 
the formation of phonon assisted satellite features as the electrons scatter off local vibra- 
tional modes. However, there are several significant differences between these systems, as 
we will see below. Electron transport through layered semiconductor structures exhibit res- 
onant tunneling features which coexist with continuum contributions. These are typically 
absent in molecular transistors. Furthermore, since semiconductor hereostructures are man- 
made, the specific resonance levels can be controlled by layer thickness and composition 
and are thus tunable. Moreover, the experimental current-voltage curves for semiconductor 
heterostructures are quite different from molecules, i.e. they show peak features rather than 



the steps characteristic for molecular systems. 



28] 



II. MODEL AND METHOD 



We wish to determine the transmission probability of electrons through potential struc- 
tures of arbitrary profile, with the possibility of exciting local Einstein phonon channels. 
The basic Hamiltonian for this problem, 

H = J2 <k)c\c k + £ V(x)c% + J2 +9 E S(x~ + K) cW, (1) 

k x xi Xi,k,k' 

describes electrons with creation and annihilation operators denoted by c* and c, and a dis- 
persion e(k) = h 2 k 2 /2m, propagating through a potential structure whose real-space profile 
is given by V(x). In addition, local Einstein phonon scatterers with creation and annihila- 
tion operators and b and energy Too are placed at impurity sites x«. The electron-phonon 
interaction is controlled by the coupling constant g, which has units of energy times length. 
The particular systems we have in mind are layered semiconductor heterostructures, such as 
GaAs/Al x Gai_ x As. For these systems, the use of a momentum independent electron-phonon 
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coupling constant is standard, and can be viewed as a reliable lowest order approach. [261] 

To find the electronic transmission probability we use the propagation matrix method, 
which is applied in the following way: for each step at position j in the potential profile, 
we construct a propagation matrix pi tep , and in between neighboring steps at a distance 
Lj apart, we construct a propagation matrix pf Tee - The elements of p{ tep depend on the 
boundary conditions of the electronic wavefunction at the potential step at position j. The 
matrix p^ ee is diagonal, and its elements depend on the phase picked up by the electron as 
it propagates through a length Lj between potential steps. The total propagation matrix is 
given by the product of the individual matrices: 



f. _ fd Ail siV 
r r step" free" 'r freer step 



(2) 



An example of the propagation matrix method is given in the appendix. For a system 
without phonon excitations, the propagation matrix is simply a 2 x 2 matrix. When the 
electron excites phonons as it penetrates the structure, the propagation matrix grows as 
(2n + 2) x (2n + 2), where n is the number of phonon channels. When several phonons are 
excited it becomes necessary to find the transmission probability of an electron as a function 
of energy numerically. The idea is to solve a system of linear equations of the form px = a, 
where x is the vector whose terms correspond to the transmission and reflection coefficients 
x = (t , r , t n , r n ) and a = (a , b , a n , b n ), where the coefficients a; and b\ depend on 
the initial conditions of the problem. In our problem there is no reflection as the electron 
exits the potential profile, so we can set the reflection coefficients r/ = for all I, therefore 
reducing the number of equations in the system by half: 
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(3) 



Initially all phonons are in the ground state, and therefore aj = 5jq. All that is left is 
to determine the tj, and we can do so by solving the system above using a Gauss- Jordan 
elimination. Once these terms are found, we can calculate the transmission probability: 
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n^ggfiM^, (4) 

where ki is the momentum of the electronic wave function in a channel with / phonons. With 
this approach it is possible to plot transmission probability versus energy. 



III. NUMERICAL RESULTS 

As a test of the validity of the propagation matrix method, we first examine two cases 
which have previously been studied in the literature. In Fig. 1(a), we show the transmission 
through repulsive (black) and attractive (red) delta potentials, which allow the excitation of 
two local vibrational modes at Tiuo = 1.0 eV and 2.0 eV. For comparison, we also show the case 
(blue) without coupling to these local Einstein modes. In their absence, the transmission 
increases monotonically with the energy of the incoming electrons. However, in the presence 
of inelastic scattering channels resonance features in the form of spikes and dips in T(E) 
occur at energies slight 



bound states 



15 



19 



.y below the local oscillator levels. They indicate the formation of 
261 ]. manifested by Fano features which for attractive potentials can 
completely suppress electron transmission right below the resonance energy (red line in 
Fig. 1(a)). These features arise from the strong feedback between inelastic and elastic 



scattering processes, and are easily missed in perturbative treatments 
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13]. The parameters 
261 ] to demonstrate 



in Fig. 1 have been chosen identical to previously published data 
full agreement of methods. 

As shown in Fig. 1(b), the phenomenon of polaron-type bound state formation persists 
for finite-width wells and barriers, [ijj] Here, the location of the Einstein scatterers are 
chosen at the center of the rectangular potential profiles. Bearing in mind experimentally 
relevant scales, we consider vibrational energies two orders of magnitude lower than in the 
Fig. 1(a), i.e. at fhu = 0.05 eV and 0.10 eV. In analogy to the case of delta potentials, 
resonances are observed at both energy levels. However, electron transmission is suppressed 
with respect to the case of delta potentials because of the finite spatial extent of the wells 
and barriers. 

Next, we turn to the case of electron transmission through more complex quantum well 
structures. Focusing on symmetric potential profiles, let us consider rectangular double 
barriers of length 0.4 nm, separation 0.6 nm, and height Vq — 1 eV. A local Einstein scatterer 




Energy (eV) Energy (eV) 

FIG. 1: (a) Electron transmission probability through a delta potential. The solid black line 
represents a repulsive potential with an electron-phonon coupling constant g = 0.117 eV ran, The 
dashed red line corresponds to transmission through an attractive delta potential with g = 0.078 
eV nm. In both cases, the electron is allowed to excite two local oscillator levels with energies hco 
= 1.0 eV and 2.0 eV. The blue line is for g = 0. (b) Electron transmission probability through 
a finite-width potential barrier/well of width L = 0.1 nm. The solid black line represents a 
repulsive rectangular potential of strength V$ = 2.0 eV with g = 0.117 eV nm. The dashed red 
line corresponds to transmission through an attractive rectangular potential Vo = —1.0 eV with 
g = 0.078 eV nm. The local oscillator levels are chosen to be at ?ilo= 0.03 eV and 0.06 eV, 
positioned at the centers of the rectangular potentials. 

is placed at the center of the well, as illustrated in Fig. 2(a). The vibrational energies are 
hu n = nhu with hut = 0.01 eV and n = 1, 2, 3, .... In the absence of phonon scattering, shown 
in Fig. 2(b), one observes a bound state at E = 0.358 eV which allows resonant tunneling 
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FIG. 2: (a) Symmetric double rectangular potential with barrier width 0.4 nm, well width 0.6 
nm, and barrier energy Vq = 1.0 eV. A local phonon scatterer is located at the center of the 
double barrier, (b-e) Electron transmission probabilities for various numbers of accessible phonon 
channels, (b) corresponds to the case without electron-phonon interaction (g = 0), (c) represents 
the transmission coefficient for one phonon channel, (d) for 2 phonon channels, and (e) for 40 
phonon channels. In (c-e) the electron-phonon coupling is set to g = 0.04 eV nm, and the phonon 
frequency is Two = 0.01 eV. With increasing number of phonon channels one observes the formation 
of a band. 

with unit transmission. In the following, we examine the effects of inelastic scattering on 
this resonant feature. In the presence of a phonon scatterer with a single available inelastic 
channel at hu\ = 0.01 eV (Fig. 2(c)), the bound state is split into two satellites, separated by 
approximately equal energy gaps with respect to the energy of the original bound state. Such 



7 



"side bands" have been the focus of numerous earlier studies. [5, 14, 15, 16, 17, 3, 3 
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Note that for the semiconductor double barrier structures studied here the magnitudes of 
the energy splits between these phonon assisted satellite peaks are considerably larger than 
the weak coupling result, E Q ± hu>i, where E is the energy of the resonance in the absence 
of inelastic scattering. This is due to strong renormalization of the bare electron-phonon 
coupling constant by the confinement of the electron wave function to the small well region, 
which will be discussed in more detail later on. 

Here, we wish to examine how such phonon assisted satellite features merge into an 
impurity band with increasing number of available inelastic channels. The generalized prop- 
agation matrix method is particularly suited for this task, as the propagation matrix for 
the system only increases linearly with the number of added vibrational modes. The pat- 
tern which emerges from Fig. 2 is that the bound state splits into n + 1 peaks, where n is 
the number of phonon channels which are excited. For instance, in the case of one excited 
phonon with energy 0.01 eV and electron-phonon coupling g = 0.04eV • nm, we find the 
peaks to be at positions E\ = 0.2833 eV and E 2 = 0.4340 eV (Fig. 2(c)), which differ from 
the zero-phonon case (Fig. 2(b)) by Ai = —0.0747 eV and A2 = 0.076 eV. For two phonons 
(Fig. 2(d)), one finds 3 peaks at E 1 = 0.2255 eV, E 2 = 0.3645 eV, and E 3 = 0.4875 eV, giv- 
ing shifts of Ai = -0.1325 eV, A 2 = 0.0065 eV and A 3 = 0.1295 eV. This observation points 
to an interesting even-odd effect, whereby for odd numbers of phonon channels, there exists 
a central, non-bonding, peak, whereas for even numbers of phonon channels it is absent. It 
also implies that the satellites at E n ~ Eq ± Hu n are bonding/ anti-bonding pairs. Before 
investigating this aspect of multi-phonon-assisted resonant tunneling more closely, let us 
point out that in the limit of many phonon channels (Fig. 2(e)) a low-energy band emerges. 
Note that the asymmetry in this impurity band is already anticipated in the asymmetry of 
the satellites for the few-phonon cases 

In Fig. 3 we examine the effects of the electron-phonon coupling and the vibrational 
energies on resonant transmission through the same double barrier potential shown in Fig. 
2(a), i.e. we focus on the low-energy transmission peaks. The case of one available phonon 
channel is studied in Figs. 3(a) and (b), and the case of two phonon channels is illustrated in 
Figs. 3(c) and (d). Let us first keep the electron-phonon coupling constant at g = 0.04eV-nm, 
and vary the energy of the vibrational levels. As observed in Figs. 3(a) and (c), increasing 
values of fkv cause the entire spectrum of transmission resonances to shift to higher energies, 
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FIG. 3: Electron transmission resonances for a double rectangular potential barrier with a phonon 
scatterer located at the center. The parameters are chosen identical to Fig. 2, unless otherwise 
specified, (a) and (b) correspond to a system with one phonon channel, and (c), (d) to a system 
with two phonon channels. In (a) and (c) the electron-phonon scattering strength is kept constant 
at g = 0.04 eV nm, and the phonon energy two is varied. One observes that the entire spectrum 
shifts to larger values of energy as hco increases. In (b) and (d) hco = 0.01 eV and g is varied. In 
this case the gaps between the transmission resonance peaks widen with with increasing g. For 
two phonon channels, the central peak does not shift as g is varied. 

whereas the gaps between the peaks remain constant. If instead we keep the vibrational 
energies fixed, i.e. at Tiu = 0.01 eV, and vary g, the gaps between peaks are found to 
increase as g increases (see Table 1). Note that for the case of even numbers of phonon 
channels the central non-bonding peak does not shift with increasing g, whereas as the 
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bonding/antibonding peaks move to higher and lower energies respectively. 



TABLE I: Gap between bonding and antibonding peaks for different values of electron-phonon 
coupling g 



g(eV • nm) 


0.02 


0.04 


0.06 


A (1 phonon) 
A (2 phonons) 


0.074 eV 
0.104 eV 


0.147 eV 
0.209 eV 


0.221 eV 
0.314 eV 



The observation of bound state energy splitting when there are Einstein phonon channels 
in the system is analogous to degeneracy breaking in the linear Stark effect. For sufficiently 
small electron-phonon coupling, we can compute the first order energy shift quantitatively 
by treating the phonon energy and electron-phonon interaction terms in the Hamiltonian (1) 
as perturbations and by using time-independent degenerate perturbation theory to calculate 
the resulting energy shifts. The unperturbed eigenstates are denoted by \x,n), where n is 
the phonon quantum number and x denotes the electron position in the well region of the 
potential profile x G [0, L]. In order to make analytical progress, the electron wave function 
is approximated by the infinite well wavefunction ty(x) ~ sm(jrx/L)/^~L, where L is the 
length of the well. The resulting perturbation matrix has elements {x,n\^ x .Twjb x .b Xi + 
g Y^ Xi ,k,k' $( x ~ x i) (&xi + ^xi) c I c fc'l x 5 n )> "which for the case of one phonon channel yields the 
2x2 perturbation matrix 

p= ( g/L 
K g/L hid 

Assuming that the impurity site is located at xq = L/2, we have sm(nx /L) = 1. The 
energy shifts are calculated by diagonalization of the perturbation matrix and are given by 



(5) 



2 



(6) 



In practice, the well width L can be made rather small, even compared to the scale of 
molecular junctions. This can lead to a significant renormalization of the electron-phonon 
coupling constant, g — > g/L, which in turn explains the the relatively large energy gaps 
between the phonon assisted satellites, observed in the propagation matrix results. 

To further illustrate how the bound state energies depend on the phonon energy and the 
electron-phonon coupling, we plot the lower and higher bound state energies as a function 



10 



of Too (Figs. 4(a) and 4(c)), and as a function of g (Figs. 4(b) and 4(d)). 
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FIG. 4: Bound state energies in a double rectangular potential barrier with one phonon scatterer at 
the center (as depicted in Fig. 2(a)). For the case of one phonon channel, the energies of the two 
resulting transmission resonance peaks are plotted as functions of the electron-phonon coupling 
constant g and of the phonon energy Ku>. Exact numerical results (solid black lines) are compared 
with the effective theory (dashed red lines) described in the text. In (a) and (c) the electron- 
phonon coupling is kept constant at g = 0.04 eV nm, and in (b) and (d) the phonon energy is kept 
constant at fiuj = 0.01 eV. In (a) and (b), for the higher bound state energies, the effective theory 
reproduces remarkably well the numerical results. In (c), the effective theory and numerical results 
are off-set by approximately 0.02 eV. In (d), the effective theory matches the numerical results for 
small values of g. For instance, when g = 0.04 eV nm, the effective theory overestimates the lower 
bound state energy by 0.01 eV. 



The accuracy of the effective theory compared to the numerical results of the full prop- 
agation matrix calculation is striking, in particular for predicting the higher bound state 
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energy (Figs. 4(a) and 4(b)). For the lower bound state, the accuracy increases for smaller 
g (Fig. 4(d)), although for a fixed value of g = 0.04eV ■ nm and variable tkv the effective 
theory predicts a lower bound state off-set by about 0.02 eV with respect to the full prop- 
agation matrix calculation (Fig. 4(c)). This difference is of the order of the phonon energy 
and one order of magnitude lower than the bound state energy. However we notice that 
although the curves for the effective theory and numerical results are off, they present the 
same qualitative behavior. Therefore, one can affirm that the effective theory reproduces 
the numerical results with striking accuracy for small values of g. The same procedure can 
be repeated for any number of phonon channels with similar results. 

Finally, let us turn to the current-voltage characteristics caused by phonon assisted trans- 
mission features. As shown in Fig. 5(a), application of an external electric field yields a 
spatial gradient in the potential energy profile. The resulting current flow is determined 
from the transmission functions T(E, V) at given voltage biases V via an integral 



where the energy window [0, Vo] for currents through semiconductor heterostructures is small 
compared to molecular junctions. As a result of this, the I(V) dependence shown in Fig. 
5(c) inherits the peak structure of the individual transmission curves, some of which are 
shown in Fig. 5(b). This is an important difference from the step-like I(V) curves reported 
in measurements of molecular junctions, s), [^j 



IV. CONCLUSION 

In summary, we have investigated how quantum well electron-phonon resonances are 
affected by the presence of several inelastic channels in the phonon spectrum. Using a gen- 
eralized propagation matrix method for multiple eleastic and inelastic transmission channels, 
we determined the highly non-perturbative effects of scattering by Einstein phonons on the 
electron transmission through potential structures. In particular, we observed a characteris- 
tic splitting of the bound state resonances into satellite peaks. The presence or absence of a 
non-bonding resonance reflects the parity associated with even vs. odd numbers of accessible 
inelastic channels. Furthermore, in the limit of many available channels, the formation of 
low-energy impurity bands was observed. The dependence of the resonance satellites on the 




(7) 
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FIG. 5: (a) Potential double barrier with thickness of 1.0 nm and separation of 5.0 nm. The height 
of the barrier is 1.0 eV, and the electronic effective mass is 0.07m e . (b) Transmission vs. energy 
curves for one excited phonon at the center of the potential well. In the presence of a voltage bias 
across the heterostructure, this curve is shifted towards lower energies, (c) Current-voltage curve 
for the double barrier in (a). The two small peaks at low voltage bias (V;,) correspond to the two 
low energy peaks in (b). The calculation of current is done by integrating the transmission over 
an energy window from eV to 50 meV. 

electron-phonon coupling strength and on the phonon energies could be reproduced using 
an effective model, which works well within the limits of perturbation theory. One promise 
of the multi-channel propagation matrix method which is developed here lies in the ability 
to study highly asymmetric quantum systems with strongly interacting itinerant and local 
features. A further direction to pursue is to depart from strictly local oscillators, which 
are nevertheless important for nanoelectronics, and to consider spatially extended phonon 
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scattering regions. 
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VI. APPENDIX 

As an illustration of how to construct the generalized propagation matrix in practice, let 
us consider the example of a rectangular potential barrier given by 



V(x) 



0, x < a 

V , a<x<b (8) 
0, x > b 



The wavefunctions are written as superpositions of plane waves, 

ip n (x < a) = a n e iknX + b n e' iknX , k n = VE - nu, (9) 
i> n (x > b) = t n e iknX , k n = VE - nu, (10) 



i>n{a < x < b) = c n e KnX + d n e KnX , n n = ^/v -E-nu, (11) 

where n represents the available phonon channels. It is assumed as an initial condition that 
the incident electrons enter the potential structure from the left, i.e. ao = 1, a„^o = 0. 

To determine the transmission coefficients for the elastic (t n=0 ) and inelastic (t n ^ ) chan- 
nels we use the propagation matrix method, matching the conditions ip n j = ip n ,j+i and 
dip n j/dx = dip n ,j+i/dx at each boundary. In the absence of phonon scattering (n — 0), the 
propagation matrix p of the system is obtained by multiplication of the step matrices (at 
x = a and x = b) and the free propagation matrix for a < x < b, 

P = PstepPfreePstepi (12) 

which are given by 

X — — 

ko j (13) 
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\ 1 j- n 



And in the interval a < x < b: 



P free 







-kqL 



(15) 



with L = b — a. 

In order to find the transmission probability, we need to solve the matrix equation pt = 
a where the coefficients of the wave functions (7) and (8) are the elements of a and t, 
respectively: 



(16) 



' Pn P12 J [ t | _ | a 

v P21 P22 ) \ J \ b 

For this system without phonon channels, the transmission probability T(E) is simply: 



l*nl 2 = 



a 



Pn 



(17) 



If instead phonon scatterer centers are present, we do not have the above condition on 
the derivative, but rather integrate Schrodinger's equation around x = (from — e to +e). 
If for instance, we add a phonon scatterer at x = a in the same potential barrier, we have 
for the step and free matrices are: 
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(18) 
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P free 





e~ KoL 

D KxL 
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(20) 





e" KlL 

Once again we determine T(E) by solving pt = a. However, now we have two transmission 
coefficients t and t±, one for each phonon channel, and the transmission probability is given 
by: 



T(E) = \t \ 2 + 



k (E) 1 11 • 



(21) 
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